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We report on thermal noise from the internal friction of dielectric coatings made from alternating 
layers of Ta20s and S1O2 deposited on fused silica substrates. We present calculations of the thermal 
noise in gravitational wave interferometers due to optical coatings, when the material properties of 
the coating are different from those of the substrate and the mechanical loss angle in the coating is 
anisotropic. The loss angle in the coatings for strains parallel to the substrate surface was determined 
from ringdown experiments. We measured the mechanical quality factor of three fused silica samples 
with coatings deposited on them. The loss angle, <j>\\ (/), of the coating material for strains parallel 
to the coated surface was found to be 4.2 ± 0.3 x 10~ 4 for coatings deposited on commercially 
polished slides and 1.0 ± 0.3 x 10 -4 for a coating deposited on a superpolished disk. Using these 
numbers, we estimate the effect of coatings on thermal noise in the initial LIGO and advanced LIGO 
interferometers. We also find that the corresponding prediction for thermal noise in the 40 m LIGO 
prototype at Caltech is consistent with the noise data. These results are complemented by results 
for a different type of coating, presented in a companion paper. 

I. INTRODUCTION 

The experimental effort to detect gravitational waves is entering an important phase. A number of interferometric 
gravitational wave observatories are being built around the world jl[ ||, |[ ||] and most should be operational in the 
next few years. Plans are already being developed to operate the next generation of interferometers; crucial research 
and development is going on now to ensure that these interferometers will have the sensitivity necessary to reach 
distances at which multiple events may be detected per year || ||, |?J . 

The sensitivity of interferometric gravitational wave observatories is limited by the fundamental noise sources 
inherent in the instrument. In advanced LIGO, thermal noise from the internal degrees of freedom of the interferometer 
test masses is expected to be the limiting noise source in the middle frequency range (~ 30 - 500 Hz). This is also the 
interferometer's most sensitive frequency band. Thus, any additional thermal noise, such as thermal noise associated 
with optical coatings, will directly reduce the number of events that advanced LIGO can detect. 

The initial LIGO interferometer uses fused silica for the interferometer test masses, the beam splitter, and other 
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optics. Fused silica has been shown to have very low internal friction |g, [| [L0[ and will therefore exhibit very low 
(off-resonance) thermal noise. This property, coupled with the fact that high quality, large, fused silica optics are 
commercially available, makes fused silica a natural choice for the initial interferometer. Sapphire, which has even 
lower internal friction plL [l^] (although higher thermoelastic loss) is currently proposed as the material from which 
to fabricate the optics for use in advanced LIGO H. In addition to lower thermal noise, sapphire offers benefits due 
to its superior thermal conductivity, which, in transporting heat from the reflective surface of the test masses, allows 
a higher power laser to be used. Fused silica is under continuing study as a fallback material for advanced LIGO 
should problems arise in the development of sapphire optics. 

In order to use the test masses as mirrors, optical coatings must be applied to the surface. To obtain high 
reflectivities, multi-layer, dielectric coatings are used. Such coatings consist of alternating layers of two dielectric 
materials with differing refractive indices. The number of layers deposited determines the reflectivity. It is possible to 
use a number of different dielectric material pairs for reflective coatings, but it has been found that coatings made with 
alternating layers of Ta20s and Si02 give the necessary reflectivity while at the same time satisfying the stringent 
limits on optical loss and birefringence required for LIGO The effect of these coatings on thermal noise is only 
now being studied. 

The simplest way to predict the thermal noise is to use the Fluctuation-Dissipation Theorem pi) . It states that 
the thermal noise power spectrum is proportional to the real part of the mechanical admittance of the test mass. 
Explicitly 

M/) = ^Jmy(/)}, (i) 

where S x is the spectral density of the thermally induced fluctuations of the test mass surface read by the inter- 
ferometer, T is the temperature, and / is the frequency of the fluctuations. The quantity Y(/) is the mechanical 
admittance of the test mass to a cyclic pressure distribution having the same form as the interferometer beam intensity 
profile For LIGO, the proposed beam profile is Gaussian. The real part of the admittance, used in the theorem, 
can be written in terms of the mechanical loss angle, rca dout , of the test mass response to the applied cyclic Gaussian 
pressure distribution. To calculate the thermal noise we must therefore obtain </> rea dout- 

The loss angle ^readout depends both on the distribution of losses in the test mass and on the shape of the deformation 
of the test mass in response to the applied pressure. If the distribution of losses in the test mass were homogeneous, the 
loss angle ^readout would be independent of the deformation of the test mass. In that case, one could obtain </> rea dout 
by measuring the loss angle associated with a resonant mode of the test mass, <fi = 1/Q, where Q is the quality factor 
of a resonant mode. However, when the distribution of mechanical losses in the test mass is not homogeneous, this 
approach does not work. 

One way of obtaining Readout would be to measure it directly. This would involve applying a cyclic Gaussian 
pressure distribution to the test mass face and measuring the phase lag of the response. But such an experiment 
presents several insuperable technical difficulties and is useful mainly as a thought experiment, in which interpretation 
of the result would be simple. 

In this paper, we give the results of another kind of experiment whose results allow us to calculate </> rea dout using 
elasticity theory. The measurement process is relatively straightforward: we compare the quality factor, Q, of vibra- 
tions of an uncoated sample of fused silica to the quality factor when a coating has been applied. In order to make 
the effect easier to measure, and to improve the accuracy of the measurements, we used thin pieces of fused silica 
rather than the relatively thick mirrors used in LIGO. Our measurements show a significant reduction of the Q due 
to mechanical loss associated with the coating. 

In choosing to make the measurements easy to carry out, we necessarily complicated the interpretation of the 
results. Scaling from the results of our measurements to the prediction of ^readout takes some work. In Section || we 
describe the relationship between the measured loss angle <j>\\ and the readout loss angle </> re adout ■ We are then able, 



in Section [II. to more succinctly explain the methodology and details of our measurement process. The results of 
the measurements are described in Section IV. The implications for LIGO are discussed Section and a program of 



future work is described in Section VI. These results are complemented by similar results on Al203/Ta20s coatings 
deposited on thick fused silica substrates. Those results are published in a companion paper p5| . 

II. THEORY 

To use the Fluctuation-Dissipation Theorem, Eq. (|l|), to predict thermal noise, we need to calculate the real part 
of the mechanical admittance of the test mass. The mechanical admittance of the test mass is defined as 

Y(f) = i2irf^f- (2) 
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where F is the (real) amplitude of a cyclic pressure distribution applied to the test mass at frequency /, and x(f) is 
the (complex) amplitude of the steady state response. Choosing the appropriate pressure distribution with which to 
excite the test mass constitutes the first step in the calculation. Levin Q has argued that in calculating the thermal 
noise read by an interferometer, the appropriate pressure distribution has the same profile as the laser beam intensity 
and should be applied to the test mass face (in the same position and orientation as the beam) . In the case of initial 
LIGO the laser beam has a Gaussian intensity distribution, and a Gaussian beam profile is also proposed for advanced 
LIGO. The corresponding cyclic pressure distribution is 

IF I -2r 2 \ 

p(r,t)=p{r,t) = — cxp —5- sin(27r/t), (3) 

where ris a point on the test mass surface, r = |rj, / is the frequency of interest, and w is the field amplitude radius of 
the laser beam. (At the radius w, the light intensity is 1/e 2 of maximum). To simplify the calculation of the response 
x(/), we make use of the fact that the beam radius is considerably smaller than the test mass radius, and approximate 
the test mass by an infinite half-space. This allows us to ignore boundary conditions everywhere except on the face 
of the test mass. For the case of homogeneous loss, Liu and Thorne (l?]] have shown that this approximation leads 
to an overestimate of the thermal noise, but that for a test mass of radius 14 cm, the error is about 30% or less for 
beam field amplitude radii w up to 6 cm. 

To calculate the real part of the admittance we follow Levin and rewrite it in the form 

Re{Y(/)} = ^U (4) 

where U (/) is the maximum elastic energy stored in the test mass as a result of the excitation, and <fi is the loss angle 
of the response. Equation (0) holds at frequencies far below the first resonance of the test mass, provided 4> <C 1, and 
is obtained as follows. Under the conditions stated 

U{f) = \\F\\x{f)\, (5) 

and the response x{f) to the excitation is 

x(f) = \x(f)\exp(-i<fr) « \x(f)\ (1 - i4>). (6) 

Substituting Eqs. (g) and (|^) into Eq. (||) and taking the real part yields Eq. (^). 

The strategy is then to calculate U(f) and <j> under the pressure distribution in Eq. (||). Calculation of the loss angle 
(j> requires some care since the loss angle is specific to the applied force distribution and to the associated deformation. 
If the material properties or intrinsic sources of loss are not isotropic and homogeneous throughout the sample, 
different deformations will exhibit different loss angles. Since interferometer test masses do have inhomegeneous loss 
due to the dielectric coating on the front surface, the calculation of thermal noise depends on obtaining the value 
of the loss angle associated with precisely the response to the pressure distribution given in Eq. (^). Throughout 
this paper we will assume that losses in the substrate are always homogeneous and isotropic and that the source of 
inhomogeneous and anisotropic loss is the coating. 

The loss angle = Readout associated with the Gaussian pressure distribution can be written as a weighted sum 
of coating and substrate losses. If the loss in the coating is homogeneous and isotropic within the coating (but still 
different from that of the substrate) we can write 

^readout — TT (^substrate ^substrate ~t~ ^coating ^coating ) (?) 

where U is the maximum elastic energy stored in the sample as a consequence of the applied pressure, E/ su bstrate is the 
portion of the energy stored in the substrate, [/coating is the portion of the energy stored in the coating, ^substrate is 
the loss angle of the substrate, and boating is the loss angle of the coating. To simplify the calculation of the energies, 
we make use of the fact that the frequencies where thermal noise dominates interferometer noise budgets are far below 
the first resonances of the test masses. Thus, the shape of the response of the test mass to a cyclic Gaussian pressure 
distribution of frequency / is well approximated by the response to an identical Gaussian pressure distribution that is 
constant in time. Thus, to a good approximation, U, J7 su bstrate, and i7 coa ting can be calculated from the deformation 
associated with the static Gaussian pressure distribution 

/ x 2F f-2r 2 \ 
p{r) = — j- exp — — . 8 

7TW \ W Z I 
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Since we are in the limit where the coating is very thin compared to the width of the pressure distribution 

U coati n g &SUd, (9) 

where SU is the energy density stored at the surface, integrated over the surface, and d is the thickness of the coating. 
Similarly, [/substrate ~ U, giving 

SU d 

^readout ^substrate ~l~ JJ ^coating. (10) 

If the loss angle of the coating is not isotropic, the second term in Eq. (|l0|) must be expanded. Since the coatings 
have a layer structure, it may not be accurate to assume that their structural loss is isotropic. To address the possible 
anisotropy of the structural loss we shall use the following model. The energy density p u of a material that is cyclically 
deformed will generally have a number of terms. We shall associate a different (structural) loss angle with each of 
these terms. For example, in cylindrical coordinates 



Pu = Prr + Prd + ■ ■ ■ (H) 

where 

Pr6 = \<3retr6 (12) 



where er^ are stresses and are strains. The associated loss angles are rr , <fi r g, etc. In this paper we will assume 
that the loss angles associated with energy stored in strains parallel to the plane of the coating are all equal. This 
assumption is motivated by the observation that many isotropic, amorphous materials, like fused silica, do not 
show significantly different quality factors for many modes even though the relative magnitude of the various terms 
in the elastic energy varies significantly between the modes Jl8| . The measurements made at Glasgow and Stanford 
Universities further strengthen this assumption |l5| — those measurements show no significant variation of the coating 
loss as the relative size of the different coating energy terms changes from mode to mode. Note that since we will 
always have traction free boundary conditions for the problems considered here, we shall always have e rz = e zr = 0. 
Thus we will have loss angles associated only with the following coating energy density components 



Pv\\ 
PuX 



-if 1 rr' 



(13) 



where are the strains and u'^ are the stresses in the coating. We define the loss angle associated with the energy 
density in parallel coating strains p v \\ , as 0y, and the loss angle associated with the density of energy in perpendicular 
coating strains, p' v ±, as cf>±. The components of the energy density in Eq. ( |l3| ) integrated over the surface of the 
(half-infinite) test mass are 



W±. = f s p' u± d 2 r. 

So that finally, to account for the anisotropic layer structure of the coating, Eq. ( |Io|) is replaced by 

SU\\ d 6U±_ d 



-'readout — ^substrate 



u r " u 



(14) 



(15) 



To obtain an expression for </> roa dout we need to calculate 8U\\ , 5U±, and U for a coated half-infinite test mass subject 
to the Gaussian pressure distribution p(r) of Eq. (^). The quantities, <5i7|| and SU± involve only the stress and strain 
in the coating. The total energy involves the stress and strain throughout the substrate 

/>oo POO 

U = 7T / dz rdr {e rr cr„ + teepee + £z Z cr zz + 2e rz a rz ), (16) 
Jo Jo 

where are the strains and er^- the stresses in the substrate. To obtain the stresses and strains in the coating and in 
the substrate we must solve the axially symmetric equations of elasticity for the coated half-infinite test mass subject 
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to the pressure distribution p(r). The general solution to these eq uations for an uncoated half-infinite test mass is 
given by Bondu et al. (with corrections by Liu and Thorne |p7| ). 

Because the coating is thin, we can, to a good approximation, ignore its presence in the solution of the elastic 
equations for the substrate. The strains in the coating should also not vary greatly as a function of depth within the 
coating, and we shall approximate them as being constant. Due to axial symmetry, e r g — £e r = e' r g = fig,. = 0. Due to 
traction-free boundary conditions, e' rz — e' zr — at the coating surface, and the same must therefore hold (to leading 
order) for the entire coating. This approximation is valid, provided the Poisson's ratio of the coating is not very 
different from that of the substrate. To obtain the non-zero stresses and strains in the coating (e' rrl e' g g, and e' zz ) we 
note that since the coating is constrained tangentially by the surface of the substrate, the coating must have the same 
tangential strains (e' rr and e' gg ) as the surface of the substrate. Also, the coating sees the same perpendicular pressure 
distribution (a' zz ) as the surface of the substrate. These conditions, which represent reasonably good approximations 
for the case of a thin coating, allow us to calculate all the coating stresses and strains in terms of the stresses and 
strains in the surface of the substrate. See Appendix |A| for the details of this calculation. 

Using the solutions for eL, cr^ , e^ , and cry derived in Appendix |A|, and substituting into Eqs. ( |l3|jT^ ), we obtain 
the required quantities 

U - (IT) 
2^/7T wY 

6U " /U = 7^ Y{1 + a')(l - a')(l - a) (18) 

sj j m 1 Y{l + a'){l-2a')-Y'a'{l + a){l-2a) 

5U±/U = ^ Y'(l-a')(l + a)(l~a) ' (19) 

(20) 

where Y and a are the Young's modulus and Poisson's ratio of the substrate, and Y' and a' are the Young's modulus 
and Poisson's ratio of the coating. Thus, from Eq. (fL5|) 



(21) 



1 d ( Y'{l + a){l-2a) 2 + Ya'{\ + a'){l~2a) 

^readout — ^substrate i , — I , t\ /-i \ *rll ' 

Vtt w \ Y {l+a')(l - a'){± — a) 

Y(l + cr')(l - 2a') - FV(1 + cr)(l - 2a) 



Y'(l-a')(l + a)(l-a) 



Substituting Eqs. ( p~T| ) and ( pl| ) into Eq. (Q) and substituting the result into the Fluctuation-Dissipation Theorem, 
Eq. (|l|), gives the power spectral density of interferometer test mass displacement thermal noise as 



o 2fc B r l-a 2 f 1 d 

^X\J ) ~~ Q/o r v 1 *rsubstrate ~r 



7T-V2/ wY Y snmc ' ^ wYY'{l-a' 2 ){l-a^) 

Y' 2 (l + a) 2 (l - 2a) 2 0„ + YY'a'{\ + a)(l + a')(l - 2a) (0„ - + Y 2 (l + a') 2 (l - 2a') ± ] J. (22) 



Equation (|22|) is valid provided that most of the loss at the coated surface occurs in the coating materials themselves 
and is not due to interfacial rubbing between the coating and the substrate, or to rubbing between the coating layers. 
If a large proportion of the loss is due to rubbing, the coating-induced thermal noise will not be proportional to the 
coating thickness as indicated in Eq. (p2[). Rather, it may be proportional to the number of layers and may be very 
dependent on the substrate preparation. 



The limits of Eq. (22) agree with previous results. In the limit that </>|| = the YY' term disappears and the 
result agrees with the result of Nakagawa who has solved the problem for that case by a different method (2(J . The 
limit of Eq. (^) in the case Y' =Y, a' = a, and </>j_ — agrees with the result obtained previously 

2fc B Tl_a2J 2 (l-2a) d \ 

S *(f) = ^JJJ (Substrate + ^= ~ 4>\\ ] • (23) 

For the case of fused silica or sapphire substrates coated with alternating layers of Ta20s and SiOg, the Poisson's 
ratio of the coating may be small enough ( < 0.25) that for likely values of the other parameters, Eq. ( p2| ) is reasonably 
approximated (within about 30%) by the result obtained by setting a — a' = 

S *W = ^TJ {^ ub ™ + 7f t ( Y ^ + Y' 0± ) } ■ (24) 
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Equation (|24|) highlights the significant elements of Eq. (|22|). It shows that in order to estimate the thermal noise 
performance of a particular coating, we must know all of Y, Y', <p^, and <f>j_. It also shows that if tf>u rs cf>±, then 
the lowest coating-induced thermal noise occurs when the Young's modulus of the coating is matched to that of the 
substrate. If Y' 76 Y, one of 4>\\ or <j>±_ will be emphasized and the other de-emphasized. This is particularly worrisome 
for coatings on sapphire substrates, whose high Young's mod ulus means that for most coatings <f>±_ is likely to be 



the main contributor to the coating thermal noise. Section III describes ringdown experiments on coated samples in 
order to determine </>ii . Unfortunately we don't obtain <f>±_ from ringdown experiments of samples with coatings on 
the surface. Since the coatings experience free boundary conditions they are not greatly compressed perpendicular to 
the surface (there will be some small amount of compression due to Poisson ratio effects). Therefore 4>± cannot be 
easily measured in such experiments, and no measurement of <j)±_ exists at the present time. Because of this, we can 
only obtain very rough estimates of the coating-induced thermal noise. We will set 4>± = </>|| > but the accuracy of our 
thermal noise estimates will remain unknown until d> 1 is measured. 



III. METHOD 



In order to estimate the coating loss component <f>u , we made measurements of the loss angles of fused silica samples 
with and without the Ta20s/Si02 high reflective coating used in LIGO. A standard way of determining the loss angle 
at the frequency of a particular resonant mode is to measure its ringdown time, r„. This allows the calculation of the 
mode's quality factor Q, through 

Q = 7r/ n T n , (25) 

where /„ is the frequency of the resonant mode. The loss angle at the resonance frequency is the inverse of the mode's 
quality factor 

4>{f n ) = l/Q. (26) 

Because of the free boundary conditions no energy is stored in strains having perpendicular components. The loss 
angle C oatcd of a resonating sample after coating is therefore related to the loss angle </> U ncoatcd of the same sample 
before coating by 

5U\\ d 

^coated = ^uncoatcd H <fi\\ (27) 

where U is the energy stored in the resonance. Similarly, as in Section O, the quantity SUn is the resonance energy 
stored in strains having no component perpendicular to the surface 

SU\\ = [ d 2 r ]T P ' lv (28) 
where S is the coated surface of the sample, z is the direction perpendicular to the surface, and 



Just as in Section [Tl], <f>\\ in Eq. ( p7j ) is the loss angle associated with energy stored in strains in the plane of the 
coating. Because we assume that all in-plane loss angles are identical, the loss angle <fi\\ is the same <j>\\ as in Section ||, 
and once measured, can be substituted directly into Eq. (p2[). 

For each sample resonance that was found, boated and </>uncoatcd were measured by recording the Q with and 
without an optical coating, respectively. The quantity (SUu d/U) was then calculated either numerically or analytically, 



allowing Eq. (27) to be solved for <j>n based on the measured values of floated an d </>uncoatcd- The resulting value for 
(j>\\ was then substituted into Eq. ( |22| ) to obtain an interferometer thermal noise estimate. 

In order to reduce systematic errors in the Q measurements, we took a number of steps to reduce excess loss 
(technical sources of loss, extrinsic to the sample) |2^, ^3|. All the Q's were measured in a vacuum space pumped 
down to at least 1 x torr, and more typically 2 x 10~ 6 torr. This reduced mechanical loss from gas damping. 

During the Q measurements, the samples were hung below a monolithic silica suspension made by alternating a 
massive bob of silica with thin, compliant, silica fibers. The suspensions and samples are shown in Fig. |l|. (The 
suspension is of the same style used previously in E3|, p3| , and 0).) The piece of fused silica rod at the top of the 
suspension is held in a collet which is rigidly connected to the underside of a thick aluminum plate supported by three 
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aluminum columns. Between the piece of rod held in the collet and the sample was a single fused silica isolation bob. 
Its function was to stop vibrations from traveling between the sample and the aluminum optical table from which it 
was suspended. The size chosen for the isolation bob depended on the sample, with the heavier sample requiring a 
larger bob. The two fibers in the suspension were monolithically pulled out of the neighboring parts using a H2-O2 
torch. These fibers had a typical diameter of roughly 100-200 /im. The normal modes of the sample were excited 
using a comb capacitor | j25| . This exciter was made from two copper wires sheathed with Teflon, each having a total 
diameter of about 1/2 mm. The two wires were then wrapped around a ground plane and placed about 1 mm from 
the face of the sample. Special care was taken to ensure that the exciter and the sample did not touch at any point. 
The position of the exciter is shown in Fig |l| Alternating wires of the comb capacitor were given a 500 V DC voltage 
while the other wires were held at ground to induce a polarization in the glass sample. To reduce any eddy current 
damping p2fl and to reduce the probability that polarized dust could span the gap between the sample and the exciter, 
the exciter was always kept more than 1/2 mm away from the sample. An AC voltage at a resonance frequency of the 
sample was then added to the DC voltage to excite the corresponding mode. Once the mode had been excited ("rung 
up" ) to an amplitude where it could be seen clearly above the noise, both the AC and DC voltages were removed and 
both exciter wires were held at ground. The sample was then allowed to ring down freely. 

The amplitude of excitation in the sample was read out using a birefringence sensor j2(| |27| or (in the earliest 
measurements) by a shadow sensor. For the birefringence sensor, a linearly polarized beam is passed through the 
sample at or near a node of the resonant mode under study. Modally generated stress at the node induces birefringence 
in the glass, which couples a small amount of the light into the orthogonal polarization, phase shifted by n/2. Thus, 
the light exiting the sample is slightly elliptically polarized. The beam is then passed through a A/4 wave-plate 
aligned with the initial polarization. This brings the phases of the two orthogonal polarization components together, 
converting the elliptically polarized light to a linear polarization that is rotated slightly compared with the initial 
polarization. The rotation angle is (to first order) proportional to the modal strain, and is measured by splitting 
the beam with a polarizing beamsplitter and monitoring the relative intensity of light in the two channels. This was 
done with two identical photodiodes and a differencing current-to- voltage amplifier. The output voltage oscillates 
sinusoidally at the resonant frequency in proportion to the modally induced strain. This signal is sent to a lock-in 
amplifier to demodulate it to a lower frequency, and the data is collected on a PC. The ringdown time r n was obtained 
by fitting the acquired signal to a damped sinusoid, or by extracting the envelope of the decay. Both approaches 
yielded the same results, although the accuracy of the former was less sensitive to corruption from noise. A schematic 
drawing of the optical readout system is shown in Fig. ^. 

For the shadow sensor, an LED is used to cast the shadow of the fused silica suspension fiber onto a split photodiode. 
The LED/diode pair is positioned close to where the suspension fiber is welded to the edge of the sample. The fiber 
near the weld point will faithfully follow the motion of the edge of the sample. As the sample resonates, the fiber's 
shadow moves back and forth on the photodiode at the same frequency. The amount of light falling on each half of 
the split photodiode changes proportionally. The currents from each half of the photodiode are then compared with 
a differential current-to-voltage amplifier as in the case of the birefringence sensor. The data acquisition and analysis 
were identical for both sensors. 

For relatively rigid but transparent samples like the ones used here, the birefringence sensor is significantly more 
sensitive and much easier to use than the shadow sensor. The shadow sensor is better suited to more compliant 
samples. In both cases however, the dominant sources of broadband noise were laser noise and noise from the 
differential amplifier. 

The samples were coated by Research Electro-Optics Corporation (REO) of Boulder, Colorado, USA. The coating 
was a dielectric optical coating consisting of alternating layers of Si02 and Ta20s. The coating was laid down using 
Argon ion beam sputtering, followed by annealing at 450°C. We chose to examine this particular type of coating 
because it is the one used on the initial LIGO mirrors that are currently installed at the LIGO sites. This coating is 
also the type currently proposed for advanced LIGO optics. 

The first samples we studied were three rectangular prisms in the shape of microscope slides (7.6 cm x 2.5 cm x 
0.1 cm) made of Suprasil 2 brand fused silica from Heraeus Quartzglas GmbH of Hannau, Germany. The surface of 
these samples was treated with a commercial polish to a scratch/dig specification of 80/50. There was no specification 
on the overall flatness or the surface figure. Two of the three slides (Slide A and Slide B) were coated on both sides 
with a reflective Ta205/Si02 coating of 3% transmittance for normally incident, 1 /im wavelength light. The third 
slide (Slide C) was left uncoated as a control. Slide A was suspended from a corner, which had remained uncoated 
due to being supported at those points during coating. Therefore, welding the suspension fiber to the corner did not 
induce visible damage to the coating. Slide B on the other hand was suspended from the center of one of its short 
edges. During the weld, the coating near the suspension point was visibly damaged in a small crescent shaped of 
radius 2 mm surrounding the suspension point. This damaged region was etched off using hydrofluoric (HF) acid. 
Table | shows the modes and quality factors for which Q's were repeatably measured. A preliminary version of these 
results was reported at the Third Edoardo Amaldi Conference On Gravitational Waves JmJ . 
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After measuring the Q's of the slides, we obtained from Zygo Corporation of Middlefield Connecticut a disk of 
Dynasil brand fused silica, 164.85 mm in diameter and 19.00 mm thick. In an effort to determine the effect of 
surface preparation on the loss due to optical coatings, this sample was made with strict specification on surface 
flatness, scratch/dig, and surface roughness. The coated surface had a surface flatness less than A/20 (A = 633 nm), 
a scratch/dig of 60/40, and a surface roughness less than 4 A rms. The back surface had a surface flatness less than 
A/6, a scratch/dig of 60/40, and a surface roughness less than 4 A rms. These specifications are nearly as stringent 
as the actual requirements for LIGO mirrors. To avoid destroying the surface with welding, an "ear" of fused silica 
was bonded onto the back surface using hydroxy catalysis bonding (silicate bonding) p8| . This ear is shaped like 
a rectangular block with a pyramid on one face. One face of the block is bonded to the sample, so that the tip 
of the pyramid faces radially. This allows the monolithic suspension to be welded with a torch to the tip of the 
pyramid without heating the sample very much. See Fig. |^. Once hung, the Q of the sample was measured using the 
birefringence readout. 

Due to the thickness of this sample (required to meet the flatness specification), only one normal mode had a 
frequency below 5 kHz. The useful bandwidth of the high voltage amplifier that was used to drive the exciter is about 
5 kHz, so measurements were possible only on this mode. This was the "butterfly" mode, with two radial nodal lines 
(£ = 2) and no circumferential nodal lines (n=0) p9| . 

After measuring the Q of this uncoated sample, it was sent to REO to be coated. It received a high reflective (HR) 
coating on one side having 1 ppm transmittance and optimized for a 45° angle of incidence. The sample was then 
rehung and the Q remeasured. As can be seen from Table O, the coating caused a significant reduction in the quality 
factor. To rule out possible excess loss due to the suspension, the sample was then removed and again rehung. During 
this hanging attempt (between successful hangings numbers 3 and 4 in Table ||), the isolation bob fell and sheared 
off the bonded ear. The bond did not give; rather, material from the sample pulled out along with the ear. A second 
ear was re-bonded at 180° to the original ear. Unfortunately, this ear was also sheared off in the same way during the 
attempt to suspend the sample. This time the source of the break occurred along the bonded surface, although some 
of the substrate pulled away as well. Finally, a third attempt succeeded with an ear bonded at 90° to the original 
ear (hanging number 4 in Table Despite the broken ears, the quality factor of the coated disk did not change 
significantly. The results of all Q measurements on the disk are shown in Table O. 

Since it is difficult in any measurement of high Q's to completely eliminate the extrinsic, technical sources of loss 
(excess loss), the quality factors measured for a given sample varied slightly from mode to mode or within a single 
mode between different hangings. Since excess loss always acts to reduce the measured Q, the best indicator of the 
true internal friction of a sample is the quality factor of the highest Q mode over all modes and hangings. The spread 
of measured Q's within single hangings was relatively small. For example, the three Q's measured in hanging number 
two (sample uncoated) were all between 3.1 x 10 6 and 2.8 x 10 6 . The twelve Q's measured in hanging three (sample 
coated) were all within 1.28 x 10 6 to 1.09 x 10 6 . As can be seen from Tables | and ||, the measured Q's also did not 
vary much between modes or hangings, nor between samples in the case of the two coated slides. The reproducibility 
of the Q's of the disk argues strongly that neither the silicate-bonded ear nor the broken ears affected the loss of the 
sample. The range of measured Q's for nominally similar situations is indicative of the level of the variable excess 
loss. Thus, for all our samples, the large difference in Q between the coated and the uncoated measurements must be 
due to the coating, and not to statistical variation, excess loss, nor, in the case of the disk, to the broken ears. 



IV. RESULTS 



Using the procedure described in Sec. [II, we obtained Q values from both the slides and the thick disk. To 
calculate <p\\ from the measured Q's we need to know the value of 8U\\ d/U for each measured mode of the samples. 
For transverse bending of the slides, the strain is approximately 

( 9 2 u y (z) 







■y 



l = J=Z 

otherwise 



(30) 



where z is the coordinate in the slides' longest dimension, y is the coordinate in the slides' shortest dimension with 
the origin in the center plane of the slide, and u y (z) is the transverse displacement of the center plane of the slides 
due to the bending. Displacements in directions other than y are zero for the transverse bending modes. This gives 



SU\\ d 
U 



= 7.2 x 10" 



(31) 



slide 



for all transverse bending modes of the slides. The butterfly mode of the disk is more complex, and an analytical 
expression for strain amplitude e(r) was not found. We made an FEA model of this sample and calculated t/ coa ti n g/£^ 
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numerically. This resulted in a value of 



'SU\\ d 




^coating 


U 


disk 


u 



5.3 x 10~ 3 



(32) 



disk 



for the butterfly mode of the disk. 

The quantities needed to calculate boated from Eq. (27) are shown in Table III. Substituting t he Q measurements 
from Table | into Eq. ( |26| ) to get the loss angles, then using Eq. ([H]) and the values in Table III in Eq. ( |27| ) and 
solving for <pu , we get 



0||, s iidc = 4.2 ± 0.3 x 1(T 4 . 
Similarly, from Eq. ( |32"|) and the disk Q's in Table || we get 

0,, disk = l.O±O.3x 10- 4 . 



(33) 
(34) 



The agreement in order of magnitude between these two measured values for <f>u sets a scale for coating thermal 
noise. This allows us to make rough estimates of the effect of coating thermal noise on advanced LIGO. The value of <fi\\ 
for the polished disk agrees within its uncertainty with the value measured for coating loss by the Glasgow/Stanford 
experiment fig] , despite the use of a different coating material in that experiment (Ta205/Al203 as opposed to 
Ta205/Si02). This suggests that the substrate surface-polish, which is of similar quality on the disk and on the 
Glasgow/ Stanford samples but less good on the slides, may be an important factor contributing to the loss. Further 
tests and comparisons are necessary before definitive conclusions can be drawn on this issue. 



V. IMPLICATIONS 

Using Eq. (^4|) for the thermal noise due to the coated mirrors, we can now estimate the thermal noise spectrum 
of the advanced LIGO interferometer. We calculated the range of coating thermal noise in the pessimistic case using 
the 0|| = 4 x 1CP 4 (from the slide results) and in the more optimistic case using cj)^ = 1 x 10~ 4 (from the disk result). 
In both cases, we assumed a beam spot size of 5.5 cm, which is the maximum obtainable on fused silica when limited 
by thermal lensing effects |3^]. We have extrapolated our results to sapphire substrates using the known material 
properties of sapphire, even though we didn't measure coating loss directly on sapphire. (There have been recent 
measurements of <pii for REO coatings deposited on sapphire Those results are in rough agreement with the 

measurements described here.) As mentioned before, the thermal noise estimates will be least accurate for sapphire 
substrates because sapphire coating thermal noise is likely to be dominated by 4>± which has not been measured. The 
Young's modulus of sapphire is considerably higher than both Ta20s and Si02 in bulk, so it seems likely that the 
coating Young's modulus is considerably less than sapphire's. We are aware of a single, preliminary measurement of 
Y' for a Ta20s/Si02 coating |32| which suggests that the Young's modulus of the coating is roughly equal to that 
of fused silica. We know of no measurements of the coating's Poisson ratio. For the purposes of estimating coating 
thermal noise, we will set the Young's modulus and Poisson ratio of the coating equal to that of fused silica. 



Table IV compares the thermal noise estimates for the four cases considered (optimistic estimates and pessimistic 
estimates on both fused silica and sapphire substrates) to the thermal noise estimates when coatings are not taken 
into account. The corresponding noise spectra for advanced LIGO are shown in Figs. [I] and These were generated 
using the program BENCH 1.13 |33| and show both the total noise and the contribution from the test mass thermal 
noise. The curves for the total noise were generated using the noise models and parameters from the advanced LIGO 
systems design document The figures show that coating thermal noise is a significant source of noise in the 

frequency band ~ 30-400 Hz for fused silica test masses and ~ 40-500 Hz for sapphire test masses. 

These estimates are only preliminary indications of the level of coating induced thermal noise. The largest source 
of uncertainty in these thermal noise estimates is that no measurement has been made of </>j_. Also, the Young's 
modulus of the coating material has not been definitively measured. The half-infinite test mass approximation adds 
extra uncertainty and this estimate needs to be refined by taking the finite size of the coated test mass into account. 
In addition, there remains the possibility that the loss associated with the different terms in the energy density p\\ 
are not equal as supposed here. However, if this were the case, the apparent consistency of the loss between different 
modes of the samples measured at Glasgow and Stanford [|l5| would be spurious. 

We have also examined the effect of coating thermal noise on the expected sensitivity of the initial LIGO interfer- 
ometers that are currently being commissioned. In initial LIGO, shot noise will be greater than in advanced LIGO 
and seismic noise will be significant up to about 40 Hz. Due to the higher level of these other noise sources, test mass 
thermal noise was not expected to be a large contributor to the total noise fy] . The addition of coating thermal noise 
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raises the overall noise in the most sensitive frequency band, around 200 Hz, by only 4% . Thus, coating thermal 
noise should not significantly impact the sensitivity of initial L1GO. 

In addition to the interferometers used for gravitational wave detection, there are a number of prototype interfer- 
ometer within the gravitational wave community. We have examined data from one of these — the 40 m prototype 
located at Caltech j$4). In this interferometer , the beam spot size was 0.22 cm and the highest Q seen for a mirror 
mode was Q max = 8.1 x 10 6 |||. Using Eq. (||) with <j>\\ = 4>± = 1 X 10~ 4 and Substrate = l/Qmax in Eq. (§4|) yields 
a predicted thermal noise of ~ 2 x 10~ 19 m/VSz at 300 Hz. This is consistent with Fig. 3 of J34| , Coating thermal 
noise is therefore a possible explanation for the broadband excess noise seen between 300 Hz and 700 Hz. The effect 
of coating thermal noise is also being explored in the Glasgow 10 meter prototype, the Thermal Noise Interferometer 
(TNI) at Caltech jjfj) and in the LASTI prototype at MIT. 

VI. FUTURE WORK 

The measurements and predictions described here indicate that mechanical loss associated with dielectric optical 
coatings may be a significant source of thermal noise in advanced LIGO. Plans are underway for experiments that will 
allow us to better understand and, perhaps, reduce the coating thermal noise. A program of loss measurements on 
various optical coatings deposited on both fused silica and sapphire substrates has begun so that the most appropriate 
coating may be found. There are also plans to try and correlate the loss angle of the coating with other methods of 
interrogating its structure. To improve the coating thermal noise without major changes to the optics, the coating loss 
must be reduced. Study of different dielectric materials is clearly warranted, and changes in the deposition process 
or post-deposition annealing might also lead to improvements. An agreement has been reached between the LIGO 
laboratory and two optical coating companies to engage in such research. 

Two main models exist for understanding the source of the excess loss in the coating. One is that the internal 
friction of the coating materials, thin layers of and Si02, is high. The other model is that the excess damping 

comes from rubbing between the layers, and between the coating and the substrate. Experiments are underway to 
test these models. 

Measurement of the unknown parameters in Eq. ( p2| ) are crucial. As discussed in Sec. [n| ringdown Q measurements 
can not determine <p± due to the boundary conditions on the free vibration of a sample. A variation of the anelastic 
aftereffect experiment [^7j , which will measure the relaxation rate of the coating after being stressed perpendicularly 
to the substrate, is being pursued at Caltech |3?j]. This experiment should give a direct measurement of <p±. 

As seen in Eq. (^2|), the coating thermal noise in an interferometer is a strong function of the laser spot size. 
Increasing the size of the laser spot reduces the effect of the coating loss on the total thermal noise, so large spot 
sizes are desirable. Large spots also help decrease the effect of thermoelastic damping in sapphire mirrors 138], so 
configurations to increase the spot size are already being considered. A spot size of about 6 cm is the largest that can 
be achieved on the 25 cm diameter test masses while still keeping the power lost due to diffraction below ~ 15 ppm. 
In the case of 25 cm diameter fused silica test masses, the largest spot size that can be achieved is about 5.5 cm, 
limited by thermal lensing p0[ . Larger diameter test masses and correspondingly larger spot sizes would be one way 
to reduce the effects of the coating loss on advanced LIGO's noise. However, this would require a re-evaluation of a 
number of advanced LIGO subsystems. 
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APPENDIX A: STRESSES AND STRAINS IN THE COATING 



We obtain the stresses and strains in the coating in terms of the stresses and strains in the surface of the substrate 
by utilizing the thin coating approximation, and assuming that the coating Poisson's ratio is not very different from 
that of the substrate. Denoting strains by and stresses by cry , this can be summarized in terms of the following 
constraints. In cylindrical coordinates, 



e' rz =e rz (Al) 

o' zz = Ozz 
<j' = O rz 



where primed quantities refer to the coating and the unprimed quantities refer to the surface of the substrate. Due to 
axial symmetry e' r0 = e' z0 = <j' rB = <r' ze = 0. We use the following relations, valid for axially symmetric deformations [ljj 

<J rr = (A + 2^)e rr + \e 6 e + Xe zz 
<r eg = Ae rr + (A + 2^)e ee + Xe zz 
cr zz = Ae rr + \e eg + (A + 2y)e zz 

(7 rz 2,fJ,6 rz 

where A and fi are the Lame coefficients. In terms of Young's modulus and Poisson's ratio, the Lame coefficients are 

A = Ya/((l + a)(l-2a)), 
M = y/(2(l + <7)). 



Combining Eqs. (Al) and Eqs. (A2), we obtain the stresses and strains in the coating in terms of the stresses and 



strains in the surface of the substrate 



/ _ A — A' / , \ , X+2jj, 

t zz — X'+2fj,' \ trr "f" t00 > "f" X'+2im> tzz 

6 ' rz = 6rz (A4) 

o' rr = (X' + 2f/)e rr + X'eee + X'e' zz 

<r' ge = X'e rr + (A' + 2fi')egg + X'e' zz 

o' zz = (J zz 
i 

O rz v rz 

where A' and // are the Lame coefficients of the coating, and A and /i are the Lame coefficients of the substrate. 

We obtain the stresses and strains in the substrate e^ , er^ from the general solutions to the axially symmetric 
equations of elasticity for an infinite half-space |l7], |l!| 

Ur(r, z) = J °° [a(k) - ^ f3(k) + /3(k)kz] e' kz .h{kr) kdk 

u z (r,z) = f™[a(k) + jf J1 f3(k)+f3{k)kz}e- kz J a (kr)kdk (A5) 
Ug(r, z) = (Axial symmetry), 



where u r (r, z) is the radial deformation of the test mass, u z {r, z) is the deformation of the test mass perpendicularly 
to the face [z being positive inward), and ug(r,z) is the transverse displacement. J\{kz) and Jo(kz) are Bessel 



13 



functions of the first kind. The functions a(k) and f3(k) are determined by the boundary conditions at the front face: 
a rz (r, z = 0) = and <r zz (r, z = 0) = p(r) |]l9|| . Using the pressure distribution p(r) from Eq. (|J) gives 



a{k) = (3{k) 



F 



An/i k 



exp 



8 



7 2 2 
• k W 



Substituting Eq. (A6) into Eqs. (A5) and performing the integrals leads to 

F(u) 



u r (r, z = 0) 
u z (r, z = 0) 



4-7r(A + n)r 
F(u) (A + 2fi) 



1 — exp 



2V2tt (A + n)/xw 



exp 



_2r^ 
w 2 

'7a I Io \ 2 



(A6) 

(A7) 
(A8) 



where Jo is a modified Bessel function of the first kind. These deformations are shown, along with the pressure 
distribution p(r), in Fig. ^. The strains in the substrate are obtained from the relations 



5u r /5r 

u r /r 

5u z /5z 

(Su z /6r + 8u r /5z)/2. 



(A9) 



These strains can now be used to find the stresses in the surface of the substrate through Eqs. (A2), and then to find 
the stresses and strains in the coating through Eqs. (A4). The results for the surface of the substrate are 



1(1- e -2r7» _ 4 e -2r7* 

F_ (± ( 1 _ e -2r 2 / W *\\ 

47r(A +AI ) \r* \ L e )) 



and for the coating 



F (jLlA— „-2r 2 /w 2 \ _ 4(A+p) -2r 2 /w 2 
2tt(A+m) \r 2 V / ™ 2 



2tt(A+m) 



F f 4 p -2r 2 /w 2 \ 
2tt ^ J 



= 



(A10) 



F(2(A+ M )-A') / 4 p -2r 2 /«, 



47r(A+/t)(A'+2ft') V w 



= 



2ir(A+At)(A'+2/j 



27r(A+^i)(A 



_ | V(A'+2//) A _ e -2r7t« 2 ^ _ 4(A'(A+ M )+V(A , +,/)) ^r 2 /™ 2 ^ 

/ //(V+2//) /-i p -2r 2 /™ 2 A , 4(A , (A+p+ M ')) 2r 2 /™ 2 N. 

'+2ft') V r2 V / 11,2 / 



(All) 



= 0. 
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Equations ( |A10|) can now be used to find the energy density in the substrate an d int egrated over the half-infinite 
volume, Eq. (|16|), to give the total energy in the substrate, Eq. (17). Equations ( All ) can be substituted into the 
expression for the energy density at the surface, Eq. (|l3|), and integrated over the surface to give the expressions for 
in Eqs. © and @. 



8U\\ and 6U 
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Slide 


Coating 


Mode 


Frequency 


Q 


A 


HR 


2 


1022 Hz 


1.1 ± 0.5 x 10 5 




HR 


3 


1944 Hz 


1.6 ±0.1 x 10 5 




HR 


4 


2815 Hz 


1.6 ±0.1 x 10 5 


B 


HR 


2 


962 Hz 


1.3 ±0.1 x 10 5 


C 


none 


2 


1188 Hz 


4.0 ±0.2 x 10 6 




none 


3 


2271 Hz 


4.9 ±0.3 x 10 6 



TABLE I: Measured Q's for transverse bending modes of the three commercially polished fused silica slides. Slides A and B 
were coated while slide C was left uncoated as a control. 



Hanging Number Coating Frequency Q 

1 none 4107 Hz 3.46 ± 0.02 x 10 6 

2 none 4107 Hz 3.10 ± 0.007 x 10 6 

3 HR (45°) 4108 Hz 1.28 ± 0.02 x 10 6 
4 1 " HR (45°) 4121Hz 1.24 ± 0.001 x 10 6 

^Ear was sheared off twice before this hanging. 

TABLE II: Measured Q's for butterfly mode of the superpolished fused silica disk. In hangings 1 and 2, the disk remained 
uncoated whereas in hangings 3 and 4 the disk had been coated. 



Test mass 






QcB 


Structural thermal noise 


material 


Coating loss 




( = 1/ 0rcadout) 


at 100 Hz, 


Sapphire 


none 




200 x 10 6 


1 x lO" 24 




0y = 1 x 10" 


-4 


15 x 10 6 


3 x 10~ 24 




011 = 4 x 10" 


-4 


4 x 10 6 


5 x 10- 24 


Fused silica 


none 




30 x 10 6 


6 x 10~ 24 




011 = 1 x 10" 


-4 


19 x 10 6 


7 x 10~ 24 




011 = 4 x 10" 


-4 


9 x 10 6 


9 x 10" 24 



TABLE IV: Comparison of structural thermal noise with and without taking coatings into account. The effective quality factor 
Q e ff (equal to the reciprocal of re adout) represents the quality factor a homogeneous mirror would need to have to give the 
same structural contribution to thermal noise as the actual coated mirror. (The effect of thermoelastic damping, important for 
sapphire, is not included in Q c b). The final column shows the strain amplitude thermal noise at 100 Hz in the advanced LIGO 
interferometer resulting from structural loss in the test mass coatings and substrates. 



Sample 


Parameter 


Value 


Units 


Slide 


Coating Layers 


14 






Coating Thickness d 


2.4 


/im 


Disk 


Coating Layers 


38 






Coating Thickness d 


24.36 


/im 


Both 


Substrate Young's modulus (Y) 


7.0 x 10 10 


N/m 2 




Coating Young's modulus (V') 


7 x 10 10 


N/m 2 



TABLE III: Physical parameters of the coating and samples. These values are used to calculate the coating loss 0y from 
Eq. ([27]). The value of the coating Young's modulus is a preliminary result from SMA/Virgo in Lyon, France. 
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FIG. 1: (a) The suspended microscope slide and exciter, (b) The suspended disk and exciter. In both (a) and (b), the entire 
structure below the steel collet is fused silica. 
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FIG. 2: Details of the attachment point. The suspension fiber is welded to the top of the ear. The ear is in turn silicate 
bonded along one of its flat faces to the uncoated side of the sample. 
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FIG. 3: Layout of the birefringence sensor. 
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FIG. 4: Strain spectrum for advanced LIGO with fused silica mirrors. The solid, straight lines represent the test mass thermal 
noise; the dashed curves show the total interferometer noise. The lighter curves were generated using optimistic assumptions 
including cf>\\ = 1 X 10 -4 . The darker curves were generated using pessimistic assumptions including <j>\\ = 4 x 10 -4 . The curve 
shown with dotted lines is the advanced LIGO noise curve without coating noise as modeled in the Advanced LIGO System 
Design document. In each case, the parameters have been optimized for binary neutron star inspiral. 
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FIG. 5: Strain spectrum for advanced LIGO with sapphire mirrors. The solid, straight lines represent the test mass thermal 
noise; the dashed curves show the total interferometer noise. The lighter curves were generated using optimistic assumptions 
including <j>\\ = 1 x 1CP 4 . The darker curves were generated using pessimistic assumptions including <fi\\ = 4 X 10 -4 . The curve 
shown with dotted lines is the advanced LIGO noise curve without coating noise as modeled in the Advanced LIGO System 
Design document. In each case, the parameters have been optimized for binary neutron star inspiral. 
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FIG. 6: Top: Pressure distribution p(r) from Eq. (|8|) (with F set to unity). The pressure distribution has the same shape as 
the laser intensity. Center: The resulting response of the surface in the axial direction, u z (r). The impression is wider than 
the applied pressure. Bottom: The resulting response in the radial direction, u r (r). As expected, the surface is pulled toward 
the center. 



